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On the ground state energy for a penetrable sphere and for a dielectric ball 
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We analyse the ultraviolet divergencies in the ground state energy for a penetrable sphere and a 
dielectric ball. We argue that for massless fields subtraction of the “empty space” or the “unbounded 
medium” contribution is not enough to make the ground state energy finite whenever the heat kernel 
coefficient a 2 is not zero. It turns out that a 2 / 0 for a penetrable sphere, a general dielectric 
background and the dielectric ball. To our surprise, for more singular configurations, as in the 
presence of sharp boundaries, the heat kernel coefficients behave to some extend better than in the 
corresponding smooth cases, making, for instance, the dilute dielectric ball a well defined problem. 
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I. INTRODUCTION 

The renewed interest in the problem of the ground state energy for the electromagnetic field in the presence of a 
dielectric body is triggered by the Schwinger suggestion Q that the Casimir effect may serve as an explanation for the 
sonoluminescence (for a review, see |^). The results by different groups are controversally Hi and although there 
is a large number of papers on this topic, basic issues are still to be clarihed. As it will turn out, the main difficulty 
is a proper formulation of the problem in the sense of finding the right physical setup. 

There are two ways to calculate the Casimir energy for a given configuration. The first way consists of summing 
up the retarded Van der Waals forces between individual molecules. The second way, which is studied in the present 
paper, is to employ the technique of quantum field theory under the influence of external conditions. (For a dilute 
dielectric ball the two methods have been shown to yield the same answer and in addition they are in agreement with 
a calculation of the Casimir energy treating the dielectric ball as a perturbation 0.) This latter way is based on the 
evaluation of closed vacuum loops with suitably modified propagators. Ultraviolet divergencies are then removed by 
means of a renormalization procedure. While in the case of conducting boundaries it is sufficient to subtract the empty 
space contribution in order to make the vacuum energy finite, this is not the case for penetrable boundary conditions, 
like the dielectric ball for instance. The main aim of the present paper is to stndy the ultraviolet divergencies for the 
electromagnetic field in the latter case. 

From the early history of the Casimir energy calculations in a dielectric using the quantum field theory framework 
we mention the papers M. and the more recent publications . More extensive bibliography can be found in 

(see also the book by Bermann |Q). 

The problem can be formulated as follows. Consider a macroscopic body whose electromagnetic properties are 
described by the permittivity e{x), which may be position dependent, i.e., which shows dispersion. The function e(a;) 
enters the photon propagator as an external held respectively as a background. Ultraviolet divergencies of the ground 
state energy of the electromagnetic held become functionals of e. In principle, it should be possible to interpret the 
ultraviolet divergencies inherent to the ground state energy in a way which yields a well dehned, unique result. This 
means that a classical model for the body must be found which allows to absorb the divergent contributions (counter 
terms) into a redehnition of some constants or some meaningful renormalisation condition must be found. 

In a standard approach one uses the heat kernel expansion in order to determine the counter terms. In (3+1) 
dimensions it is the coefficient 02 which must be known for a massless held, having in mind the electromagnetic 
held for instance. By now, it is still unknown for a dielectric ball. The problem is that such a background is not 
diherentiable and the standard recursive formulas cannot be applied. On the other hand, the known results for 
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manifolds with boundary (e.g. with Dirichlet boundary conditions ) do not apply because the dielectric ball results in 
some penetrable boundary conditions . Without attempting to solve this problem in general, we calculate the relevant 
heat kernel coefficients in two sufficiently characteristic examples, first a spherical surface carrying a delta function 
potential (’delta sphere’, V{r) ^ S{r — R)) which can be reformulated in terms of matching conditions at r = R, and 
second, the dielectric ball with constant permittivity and permeability inside and outside. In both cases we discuss 
in parallel the corresponding formulas for a smooth background. 

Before actually calculating the heat kernel coefficients we discuss the renormalization scheme in general, especially 
the problem of finding a normalization condition. We argue that there is a natural condition for the massive field and 
show the difficulties for the massless case caused by the conformal anomaly. 

Because of the contradicting results mentioned above we found it useful to reconsider the quantization procedure 
for the Maxwell field in a dielectric. In our previous paper p5[ | the quantization in terms of usual vector potentials 
was considered. We obsereved no cancellation between ghosts and “non-physical” modes. Here we re-examine 
this problem in terms of the so-called dual potentials. We give a rigorous proof that for the spherically symmetric 
dielectric ghost contribution can be neglected. We also show that the both approaches give identical results for the 
heat kernel coefficient 02 in the dilute approximation. 

The paper is organized as follows. In the next section we discuss the renormalization scheme, in sections 3 and 4 we 
calculate the heat kernel coefficients for a delta sphere and for the dielectric ball, the results are discussed in section 
5. Some formulas are given in the appendix. 


II. RENORMALIZATION AND CONFORMAL ANOMALY 


Let us consider the ground state energy of a massive quantum field in the zeta-function regularisation, 

,2s ^ 


^0 = ^ Cvis - -), 


where 
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(ri) 

is the zeta function of the operator V corresponding to a given dynamical system, 
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The index (n) includes discrete and continuous parts in dependence of the problem considered. An arbitrary parameter 
/i has the dimension of a mass. The zeta function is connected with the heat kernel by means of 
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C(s) = / dt 
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The first few terms (up to n = 2 in (3-1-1) dimensions) of its asymptotic expansion for t > 0 
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yield the ultraviolet di verg encies of the ground state energy ([^. From that we define the ’divergent part’ of the 
ground state energy as ||^ 
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Now the renormalised ground state energy can be defined as 
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in the limit s —> 0. In general, the renormalisation is not unique as it is obvious already from the presence of the 
parameter fi. Therefor some normalisation condition is necessary. 

In general, it depends on the considered physical problem how to handle this non uniqueness. The simplest example 
is the Casimir effect between two separate reflecting bodies with surfaces dAii and dA42- In flat space and without 
background potential the coefficients a„ are integrals over the surface and for 9AIi ni 9 Ad 2 = 0 they are independent of 
the distance between the bodies. Therefor all divergent and non unique contributions do not depend on the distance 
and a force as a physical quantity can be defined. Note that this is sufficient for all current experiments on Casimir 
force measurements. This applies especially to the case where dA4i and dA42 are nonintersecting spheres. But it 
does not apply to the case of a single sphere when asking for the pressure on its surface , for instance. 

In the case of a massive field there is a natural normalisation condition. It is the requirement that the renormalised 
ground state energy must vanish for the mass m becoming large, 

—!■ 0 for m ^ oo . ( 8 ) 

Because the heat kernel expansion is equivalent to the asymptotic expansion of the ground state energy for large m 
(provided there is no other dependence on m than that in Eq. (^)), the definition of the divergent part through 
the heat kernel coefficients is just equivalent to (H). This condition had been used in @ for the calculation of the 
ground state energy of a massive scalar field with boundary conditions on a sphere. 

Regrettably, this normalisation condition cannot be extended to the case of a massless field. In general, for to = 0, 
all dangerous contributions are proportional to 02 . In case 02 = 0 there is no divergence and no arbitrariness at all 
and, hence, no need for a normalisation condition. But for 02 ^ 0 the limit to —*■ 0 cannot be performed in ( 0 ) 
as can be seen quite easily. For this we observe that the ground state energy © can be expanded into a power series 
for small to provided there are no zero modes (we assume all A(„) > 0 ) 
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Now, by means of (|^) and (i we subtract a contribution containing (02/1671^) login. Therefor, the behaviour of the 
renormalised ground state energy eq. (|^), for small to is 
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Consequently, there is no transition from the properly normalised ground state energy of a massive field to the 
corresponding massless case as long as 02 0 . 

A more general procedure for handling the divergencies is the following. One considers the quantum field in the 
background of some classical system. One may think of quantised matter fields in the background of the classical 
gravitational field or of quantum fluctuations around a classical solution. The minimal structure of the background is 
determined by the heat kernel coefficient a„ with n = 0 ,..., 2 which are certain functions respectively functionals of 
the background. In fact, the structure of the energy of the classical system is determined by £o‘'', eq. (||). In writing 

^•ciass _ _j_ (jaii2 + FCLi -\- A:a3/2 + ha2, (11) 

where p, a, F, k, h are the parameters of the classical energy, the subtraction in ( 0 ) can be interpreted as a redefinition 
of these parameters. They remain undetermined and, hence, cannot be predicted from the calculation of the ground 
state energy . In fact, by means of this, the heat kernel coefficients give the physical prediction which contributions 
are of quantum origin and which are of classical nature. 

Again, a remark must be added for the massless case. Consider, for example, a scalar field with boundary conditions 
on a sphere. The dependence of the eigenvalues on the radius is A(ji) = A(„)/i?^ where A(„) are numbers (depending 
on the kind of boundary conditions ). Hence the ground state energy can be written as 


2R 


Cvis-1/2), 
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where Q-p is the zeta function for to = 0 and R = 1. In case 02 7 ^ 0, one has the expansion 

C(s - 1 / 2 ) = +h + 0 (s), 
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where 02 = a2{R = 1 ) and h is some number which cannot be calculated from the asymptotic expansion. The ground 
state energy becomes 


^" = 325r(; + ‘°«™')+S' ™ 

The appearance of logi? is just a result of the conformal anomaly. This was clearly stated in 0, where the scaling 
behavior of the Casimir energy was investigated (see also ). In order to setup the renormalisation scheme one has 
to introduce the classical energy 


^class 

R 

which allows to absorb the pole term into a redefinition of the parameter h. However, the contribution h/R in £q 
which could be viewed as the genuine result of the calculation of the ground state energy (sometime it is called the 
nonlocal contribution) cannot be distinguished from the classical part and its calculation does not have a predictive 
power. Thus the only outcome from the calculation of the ground state energy is the contribution containing logi?. 
This consideration applies in the case 02 ^ 0 . No renormalization procedure and no classical system are required in 
case 02 = 0 . The ground state energy can be calculated directly and becomes 


^0 = ( 15 ) 

We remind that for the scalar field with Dirichlet or Robin boundary conditions and for the electromagnetic field 
with conductor boundary conditions , the coefficient 02 is the same except for the sign for the interior and the exterior 
regions (the reason being that the extrinsic curvature has opposite signs on both sides). Therefore, when considering 
the ground state energy in the whole space, 02 is zero and a result like ( 0 ) can be obtained. An example for this is 
just the well known Casimir effect for a conducting sphere, where ’delicate’ cancellations have been observed in early 
papers already. In contrary, when considering only the interior (or the exterior or an odd dimensional spacetime), 
one has 02 0 and only the logarithmic contribution in (^^ is the part which it is meaningful to calculate, i.e., it is 

sufhcient to calculate 02. 

From this point of view, let us consider a conducting sphere of finite thickness. Let a massles scalar field be given 
in the regions (0 < r < i?i) U (i?2 < r < 00), Ri < R2, obeying Dirichlet boundary conditions at r = i?i and r = i?2- 
Let the field be zero in between the two spherical shells having in mind an ideally conducting medium there. The 
heat kernel coefficient is 0,2 = —167r/315. In zeta functional regularization, the ground state energy for a single sphere 
is given by eq. ( 0 - From that we obtain for two concentric spheres 
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where hi ^ /12 are the ’genuine’, nonlocal contributions of the ground state energy inside respectively outside a sphere 
of unit radius. 

Now, by writing 


f _ ^1 + ^2 

i?i i?2 “ ^ V^i R2) i?2 

we can define a classical energy 

where h is some free parameter, and the renormalised energy by (|) 
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which delivers a meaningful expression for a pressure on the surface at r = i?2- 
By writing 
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f _, ^1 + ^2 

i?i i?2 “ R2J Ri 

we get a similar result on the inner surface at r = i?i. Thus, for the given configuration, one is forced to introduce 
a classical model with the energy (0 in order to obtain a renormalised Casimir energy. However, although this 
procedure might look natural, it isn’t. So in the first case, the renormalized ground state energy is independent on 
the inner radius which seems unlikely to be reasonable. 

In summary, whereas for the massive field a unique normalized Casimir energy may be defined by considering the 
m ^ 00 behaviour, for the massless case, as seen for different situations, finite renormalizations remain undetermined 
as soon as the heat kernel coefficient 02 7^ 0. 


III. PENETRABLE SPHERICAL SHELL 


In the preceeding section we discussed the role of the coefficient 02 for the ground state energy , for instance the 
well known fact that there is a cancellation between the two sides of a surface dividing the quantisation volume. Here 
we show that this is a rather special case and that generally the renormalization ambiguity remains. We start from 
the known example of a smooth background field V(x). The operator V in is 

P = -A + V(x) ( 18 ) 


and the corresponding heat kernel coefficients are 


Oi = - 
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daf V(x), 



da? V(xy. 


( 19 ) 


The coefficient 02 can vanish only in the trivial case of a vanishing background potential. Consequently, the presence 
of 02, i.e., of the conformal anomaly, is a rather general case. The question arises how this is connected with the 
cancellation of 02 in the presence of boundaries. As is known, there is no transition in the coefficients from a smooth 
background to boundary conditions . For instance, boundary conditions lead to coefficients with half integer numbers 
which are not present in (Hi)- In the following we consider a singular background potential 

V{r) = ^S{r-R) , ( 20 ) 


which can be viewed as standing in the gap between a smooth background and boundary conditions. The delta 
function potential ( 20 ) in the operator V, eq. (|l8|), can be replaced by penetrable boundary conditions , namely by 
the matching conditions 


^\r=R-0 '^|r=R+0> *^|r=il-0 ^\r=R+0 4’\r=RJ (^l) 

on the solutions of eq. (^ with V = — A. Clearly, the fomulas (^ cannot be applied because 02 would contain the 
delta function squared. 

For a > 0 this potential is repulsive and in the formal limit a —*■ 00 Dirichlet boundary conditions at r = i? are 
restored. The Casimir force between two planes carrying delta potentials had been calculated in [p^ , po[ and the similar 
problem for dissipative mirrors in pl|] . In these cases one had been interested in the distance dependent contributions 
only. 

For the calculation of the ground state energy and the heat kernel coefficients we adopt the technique developed 
in for a smooth background potential. The starting point is the expression of the ground state energy in zeta 
functional regularisation 


00 

^ 1=0 

with V = 1 + 1 / 2 . Here, fi{ik) is the dost function of the scattering problem associated with the operator V, eq. (|^), 
on the imaginary axis. It can be found by standard methods and reads 

Mik) = 1 + aR{kR)K,{kR), ( 23 ) 
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where Ii, and Ki, are the modified Bessel functions. In order to calculate the heat kernel coefficients it is sufficient to 
calculate the pole part of 5o in s = 0. For this reason we use the uniform asymptotic expansion of the dost function 
for k —>■ oo, V ^ oo. From the corresponding asymptotic expansions of the Bessel functions |23| we obtain 




(24) 


with t = 1/+ {kR/vY nnd Xi i = a/2, X 2,2 = a^/S, ^3^3 = a/lG + a^/24, X^^^ = —3a/8, X^ j = 5a/16. Higher 
terms of the expansion do not contribute to the pole. Inserting log/®'® into the representation (|2^) we denote the 
corresponding part of the ground state energy by £®®. There the integration over k can be carried out explicitely 
(e.g., using formula (31) in ^^). After that, the sum over I will be converted into two integrals by means of 
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and £®® splitted into two parts f accordingly. In the first part the integral over ly can be carried out 

explicitely and we obtain 
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In from Aii^i after some calculation a further pole contribution results 
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where all other contributions vanish for to —> 00 and hence do not contribute to the considered heat kernel coefficients 
. Comparing these expressions with (ra we read off the coefficients 


01/2 = Oj oi = —dTrai?, 


03/2 = 


27r a® 


(25) 


(oq does not depend on the potential V{r)). 

From these formulas we see that the coefficients ai/2 and oi are the same as for a smooth background. In fact, in 
( |l^ ) ai/2 = 0 and ai equals ( ^ when inserting V{r), eq. (^0|). In contrast, O3/2 in ( ^^ has no counterpart in the 
smooth case, which is not surprising as in the case of boundary conditions coefficients with half interger numbers are 
present. Finally, the coefficient 02 in (|^ ) cannot vanish except for the trivial case a = 0 and it cannot be obtained 
from 02, eq. (H), because of a squared delta function which would appear. It is interesting to note that contributions 
linear in a, which are present in intermediate expressions and cancelled. We will see below that for the 

dielectric ball there is a similar cancellation. In order to renormalize the groundstate energy one needs to introduce 
the classical part = ha 2 , see eq. (0), of the energy and there is no condition at hand to fix the free parameter 

h. 


IV. DIELECTRIC BACKGROUND 

There is a general interest in a careful consideration of the quantisation procedure in a dielectric background. The 
ground state energy is the sum over all states of the considered Hilbert space, physical as well as unphysical ones. 
This is in opposite to the usual situation where one has to consider expectation values in states with particles, i.e., in 
states of the physical subspace. Therefor the unphysical states, including ghost states, must be considered carefully. A 
general investigation whether the contributions of the ghosts and of the unphysical photons cancel each other is still to 
be done. So a reconsideration of the in general well known quantisation procedure of QED in a dielectric background 
is in order. Thereby we consider two quantisation procedures, first in terms of the electromagnetic potentials and, 
secondly, in terms of the dual potentials. In fact, the first case was already done in [|^. We reconsider it here in 
order to derive specific expressions in the spherically symmetric case. 
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The action for the electromagnetic field in a medium with permittivity e(5’) and permeability /i = 1 is 

S=^J d*x{e{x)E^ - B^). (26) 

First we consider the case of e(x) being a sufficiently smooth function so that the necessary number of derivatives 
exists. Furthermore, we assume that e(x) goes sufficiently fast to a constant value at spatial infinity such that integrals 
appearing later on are well defined. Volume divergencies still present have to be subtracted by hand. These terms 
are associated with the contributions of an unbounded medium, where e everywhere takes its value at infinity. As is 
well known, for the calculation of the ground state energy in gauge field theories care must be taken in order not to 
loose modes, ghost contributions for instance. This had been done in |^. There, after introducing standard vector 
potentials, a gauge fixing term, ghost helds, and taking into account the proper functional measure, the path integral 
takes the form 


Z = y DAiDAQ Dc J d'^x[^{2€{x){doe ^^‘^Ai — diAo)'^ 

— {die Ak — dke + Cgf + /lg?iost]|’ 


(27) 


with Ai = \feAi and 

^gS = -\^{e~^die^^'^Ai - e^oAo)^, 

^ghost — ^ dicdi “h 6(9 q)c. 

The action for the electromagnetic field A is then found to be 

i J d'^x[e{diAof - e‘^{doAof + {doAi)'^ 

+Aie ^^^{d'^Sik — Cidk + diCk — eiek)e 

with the notation = (9dne, and later also Cij = didjlne and so on. 

Now we are able to integrate over Aq, A and the ghosts. The resulting path integral reads after Wick rotation to 
the Euclidean domain 


Z = Z[Ao]Z[A]Z[c,cI 


where the separate contributions are of the form: 

Z[Ao] = det~^^'^{-d,ed, - e^Sp), 

Z[A] = det“^^^ {^^dldij - d^dij - Ciedj + ejed^ - ^{cij - e^ej)^ , 
Z[c, c] = det(—e“^i9ie9i — ed^). 


For the functional determinants we use the integral representation 

/»00 T I 

logdet(L) = / —K{L;t), 

Jo * 

where the heat kernel K{L;t) for a second order elliptic operator L is 

K{L\t) = Trexp(—tL). 

In |l^ we calculated the coefficients an of the asymptotic expansion of the heat kernel 


Tr(exp(—fL)) = (Ant) ^ ^ Ean{L). 
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(28) 
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Note, that here L is a four-dimensional operator. In eq. (^) above 7^ is a three-dimensional operator. This explains 
different powers of 47rt in (^) and 

Let us list the individual contributions to the different coefficients (summation over repeated indices is assumed). 
First 
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For oi we obtain 
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and for the coefficient 02 the result is 
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(in the first line of we corrected some tipos which appeared in |p^ ). Because the background, which we consider 
here, is static we can drop the integration over xq. The contributions must be summed up according to 


an = a^o +a^- 2af. 


(29) 


The expression for 02 appearing in this way represents the necessary counterterm which is required to perform the 
renormalisation of the ground state energy of the electromagnetic field. As it is seen, this is a quite complicated 
functional of e(x). It is interesting to consider the case of e close to unity (dilute dielectric medium). By means of 
e = 1 + S, Bi = Si/(1 -I- (5), etc. we obtain 

/ 19 

dx (2-f35-f — 

ai = I +0{S^) 

and 

02 = ydx ^ {Suf + O(S^). (30) 


In oi and 02 there are no terms of the first order, they cannot be present for symmetry and dimensional reasons. 
(In some detail, the coefficient oi has the dimension of length, 02 of 1/length and there is no possibility to have 





nonvanishing terms linear in <5 which do not integrate away.) But the terms of the second order are present. Moreover, 
they cannot vanish (unless e is globally constant). 

From this one can conclude, that for a smooth e{x) the ultraviolet divergencies in the ground state energy of the 
electromagnetic field are present generically and can cancel at most for very specific choices of e but not for small 
e — 1 as shown by (^. 

The functional dependence of the needed counterterms is very involved and the classical model necessary for 
renormalization lacks physical intuition as well as motivation. 

So we obtained the relevant heat kernel coefficients for a e{x) which must be at least for times differentiable in order 
to deliver a finite result when inserted into the equation for 02. When restricting to a spherically symmetric situation, 
the quantization approach using dual potentials will provide several advantages. (For the general formalism of dual 
potentials see |^^.) Probably most important is, that the calculation of the relevant determinant can be reduced to 
the problem of two scalar fields. This procedure can be taken over to the non-smooth dielectric ball treated below 
and for that reason is described in some detail. 

In this approach one starts with the action 

S =-^ J d^x (31) 

which generates the Maxwell equations in dielectric medium with permittivity e{x). Here D = eE. Introducing the 
dual potentials D — V x C, B = dtC — VCo, one gets the canonical Poisson brackets, 

{Ci{x,t),Bmj,t)} = Sl6{x-y), (32) 

thus B is the momentum conjugate to C. 

One can rewrite the action (§) in the canonical first-order form, 

5 = j d^x (^BdtC + Co(VH) - - 1^2^ . (33) 

Co plays the role of the Lagrange multiplier. 

The theory (^^ is then quantized according to the general method j^] . The path integral measure reads 

dfi = DCjDB’^DCoJppdix) , (34) 

where x is a gauge fixing condition, Jpp is the Faddeev-Popov determinant, Jfp = det{x, (VH)}. For x we choose 
the condition 


X = V.C* = 0 


(35) 


The Faddeev-Popov determinant becomes 


Jfp = det(-V*V,). 


(36) 


The ghost operator in (36) is non-elliptic, and, therefore, the determinant ( |36[ ) is ill-defined. Such operators must 
be treated with great care in order not to break the gauge invariance of the theory (see e.g. ^^). Fortunately, in 
our case Jpp does not depend on e and can be neglected alltogether. Moreover, at the expence of some technical 
complications we can obtain an equivalent result in the Lorentz type gauge = 0 where no such problems occur. 

Now one can perform the intergrations over all fields arriving at the following expression for the path integral. 


Z[e] = J diiexp{iS) = j Z?C_l exp(iS'), 


(37) 


where the subscript T denotes transversal 3-vectors, iC\ = 0. 

Two substantial simplifications occur here compared to the quantization in terms of non-dual vector potentials 
described before |^. First of all, the gauge condition on Ai (^), which removes the mixing between Aq and Ai, 
depends on e. As a consequence, the ghost operator depends on e as well and cannot be neglected. Also, in the case 
of non-dual potentials the term qudratic in momenta is multiplied by e~^. After momentum integration it leads to a 
modification of the path integral measure. 

Even though the path integral ( p7] ) is Gaussian its direct evaluation is not easy. The corresponding determinant 
is to be computed on the space of transversal vectors. To use standard heat kernel methods one must extend the 
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functional space to all unconstrained vector fields at the expence of introducing a compensating scalar determinant. 
This procedure is equivalent to introducingghost fields. However, if the function e(x) is spherically symmetric, one 
can reduce the determinant arrising from (^) to two scalar ones by making use of standard TE and TM modes. We 
turn to this case, now. 

The transversality equation (^) can be solved by introducing two independent scalar modes (j) and tjj, 

Cr = rAtjj, 

Ca = -VadrTip + readb4>- (38) 

Here the coordinates x“ denote angular variables, A and V are the Laplace operator and covariant derivative on the 
two-sphere. Sab is the two-dimensional Levi-Civita tensor. 

It is clear from that the scalar harmonics with zero angular momentum do not give rise to any transversal 
vector fields. Such harmonics must be excluded from the path intergral. 

The normalization conditions for the scalar modes one can read off from the integral 


= f r^dr + (j)[—Kr^)(j^ 


(39) 


where dO is the measure of the angular integration. 
By substituting ( p^ in the action (ED we obtain 


S = / r^drdn 


i/)(r^A)A 9o + ^ 


1 


+(j){-r^A) -9q -|- 9r + - - 9r.-l-- -l-A 


(40) 

We can now evaluate the path integral (^). Taking into account the Jacobian factors appearing due to (^^, one 
arrives at the following expression 


Z[e]=det' '/^(Lv-)det' (41) 

with 

Lip = H—A, 

L(j) = —dg H—A + ^dr—^ ^dr H—^ . (42) 

The prime in (1) reminds us to subtract zero angular momentum fields. The operator is not hermitian. However, 
its heat kernel coincides with that of the hermitian operator 

The heat kernel coefficients for the operators and which include zero angular momentum fields can be easily 
calculated by using the standard methods . After continuation to the Euclidean domain the oparators (^) become 
elliptic operators of the Laplace type. This means that they can be represented in the form: 


L = - J- J- h) 


(43) 


where is a symmetric tensor field which plays the role of the metric. By introducing a connection one can 
bring L to the form: 


L =^ + E) (44) 

where V is a sum of the Riemannian covariant derivative with respect to the metric g and the connection uj. The 
explicit form of u) and E is 

cos = \g.5{a’^+ 

E = b-g^''{d^u)„ + uj^LOy - (45) 

As usual, r denotes the Christoffel connection. We must stress here that the metric g^j^y is an auxilliary object. It 
does not reflect any real geometry of the problem. The heat kernel coeffcients can be expressed in terms of the fields 
introduced above IM: 
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ao= (Txg 


1/2 




ai =y d'^xg^/'^[E + 

a2= f d'^xg^/^ — ieOrE + 180E^ + 30n^^n>^'' + - 2p^ + 2R^) 

J 360 


( 46 ) 


Here i?, p and r are Riemann tensor, Ricci tensor and scalar curvature of the metric g respectively. Semicolon denotes 
covariant differentiation, E-^^ = ^fiE. All indices are lowered and raised with the metric tensor, fl is the field strength 
of the connection w, which is zero for the operators (42). This is the most effective, though maybe a bit artificial way 
to calculate the leading terms of the heat kernel expansion on a smooth background. 

Metric and curvature tensors coincide for the operators and L^. In the dilute approximation they are: 

9ij ~ I 

E jkl — G-ildkj T . . . 

Pjk — "^i^jk 4” ^ppdjk) 4“ . ■. 


T = —2e 


pp 


where we dropped the terms which do not contribute to the first non-trivial order in the (e — 1) expansion. The 
functions E are different for and L^: 


Eib — -rGii 4- . . . , Ea, — —err 4- V—h . .. . 

4 4 2r 

From these expressions we obtain the conformal anomaly in the dilute approximation: 


(47) 


a2'[4>\ = J r^drdn^ 
02 'M = J 


39 51 129 ^2 1 

- - - 5“^ 

2 r2 4 ’ 



(48) 

(49) 


Note, that the coefficient (^) can not be represented in the covariant form of the coefficient (^^. The prime means 
that the zero angular momentum modes (0° and are still to be subtracted. Introduce the fields 0° = r0° and 
•0° = which are square integrable on the half-line with the unit weight. Acting on the fields 0° and (jr the 
operators and become 


= -dl E -dl 
L\ = -dl + drUr 


(50) 


Since the operators ( pOD are two-dimensional, the coefficient before the zeroth power of proper time is given by oi. 
Again, it can be calculated by the same methods of The result is 




1 


dr ( -.0.^ 1 = - 


1 


r‘drdnt-Ed 


(51) 


The total coefficient 02 after subtracting the zero angular momentum modes reads 

02 = a2'[0] 4- a2'[0] - 47r(a?[0] 4- a?[0]) 

= [ dr dVt^-5li , 

/ 720 " 


(52) 


where the factor 47r appeared due to different normalization conventions for the heat kernel coefficients in 2 and 4 
dimensions. We see, that the covariance of the result is restored. The coefficient (|^) coincides with the result of the 
calculations in ordinary potentials, eq. (^^, providing a check of our previous calculation. 
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Now, because the problem is reduced to two scalar ones we can turn to the case of a dielectric ball of radius R 
whose permittivity and permeability are constant inside and outside 

e(r) = ei 0(i? - r) + e2 0(r - R), (53) 

(and similar for We remind that it is just the step contained in this e(r) why the above approach cannot be 

applied. 

The most convenient form of the operators (^ ) to analyse the present case is 

L^ = -d^ + ~d^r+jA, (54) 

= —Oq + -dr-drT + - A. (55) 

r e e 

The corresponding eigenfunctions, which can be divided into two parts according to 

'ip(r) = 'ipi{r)Q{R — r) + i/)2(r)0(r — R) (56) 

(similarly for (jj), are determined by ipir) and its first derivative being continuous, 

r^|_i_ = rip\- i9rC'0|+ = drri/^l- (57) 

and for (j) we have 

r(j)\+ = r(j)\- ^drr(l)\+ = , (58) 

i.e., the function and a known combination of the first derivative with e must be continuous. It is clear that these condi¬ 
tions are just the same as known from the standard approach, see for example a textbook on classical electrodynamics 
as |§. 

Note that this situation is quite similar to the delta shell (^l|). Just like in that case, to proceed further one has to 
define the corresponding scattering problem and the corresponding Jost functions. They turn out to be proportional 
to the expressions known for example from |^] 

Af^{kR) = y'eifi2s'i{kiR)ei{k2R) - y/e^si{kiR)e'i{k2R) 

Af^{kR) = y/e^s'i{kiR)ei{k2R) - ./e^si{kiR)e'i{k2R) (59) 

with the notations 


si{x) = 




(60) 


{u = 1+1/2). Here, ei.2 and /ri.2 are the permittivity and permeability (we keep it in these formulas) inside respectively 
outside the ball. Accordingly we use the notations fci^2 = These expressions for and A^^ can be 

inserted into formula (|^ ) for the ground state energy whereby any multiplicative factors (even dependent on k as 
long as they are analytic near the real axis) by which A^^ and A^^ may differ from the Jost functions, drop out. 
By means of Eq. (i the corresponding zeta function then reads 

C(s) = ^(2^ + 1) dk — In [Af^{kR)Af^{kR)] . (61) 

As will be seen clearly afterwards (eq. dzl)). the Casimir energy of the dielectric filling all space, and having permit¬ 
tivity £2 and permeability /i2, has already been subtracted. This is the translational invariant so called ’Minkowski 
space’ contribution. 

It is useful to introduce the notation 


Ap_/(fca) = i''s'i{kia)ei{k2a) - s/(/cia)e((fc2a) (62) 

with ^ and with the correspondence Aj^ Ai^; and Aj^ A_i.i. In passing from Eqs. (^) to d^ ) 

only irrelevant factors have been neglected. 
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As a result, after substituting kR = zv^ the zeta function of TE respectively TM modes read 


Cp(s) = ^E( 2 ^ + l) / 


ZI/N- 2 S d 


l=l 




dz 


In [^p,i{^^)\ ■ 


( 63 ) 


In order to calculate the heat kernel coefficients we need the uniform asymptotic expansion of It can be obtained 
just as in the preceeding section by using the known uniform asymptotic expansions of the Bessel functions. But 
because the resulting expressions are more involved now, we introduce some obvious abbreviations involving the Debye 
polynomials Mfe(t) and Vk{t) | |^ , 

^ Ufc(t) o v^oo ! 

^ — 2 ^k=l jyfc ’ £) — lyk ) 

^ V^CX 3 Vkit) J~, v^oo ! T\kVkit) 

In terms of these the asymptotic expansion of s/, ej, and their derivatives read 


Sl{vz) = 


ei{uz) = 
s'i{vz) = 


e[{vz) = 


A /2 


2(1+22)1/4 




A /2 


(1+22)1/4 

1 1 


e-‘^''(^)(l + B), 

+A) 


Av 21 / 2(1 + 22 ) 1/4 


2+2 

1 


21 ^ 21 / 2(1 + 22 ) 1/4 

(1 + 2^)^/"^ 

zl/2 


e-^’'(^)(l + B) 


(1 + D). 


Before stating the asymptotic expansion of Ap^; some further notation is necessary. First, the velocities of light are 
defined by Ci = Ijy/eijll (i = 1,2). Furthermore we use ti = 1 /-^/l + 2^, Zi = ^ti^iZ, and we will write Ai, Bi,Ci, Di 
indicating to which variable the contribution belongs. Then, after some calculation, we get 


2 ^ Cl€2/1/2 


1 + 


V Cl €2/1/2 
+ C2tl 




(Ci + ^(1+Ai))(1+B2)+B2 


C2/1 

Cl/2 


/2 


(D2-;^(1 + B2))(1+Ai)+Ai 
2v 


The expansion for v ^ 00 oi the logarithm entering 

InAp^i - 

with 


reads 


E 

n=-l,0,l, 




D_i,p(2)=r? - -ly - 


Cl 


C 2 


Do^piz) = In - 


1 'f^Ci/2 + C2/l 


2 Cl 02/1/2 

the consecutive expressions are listed in the appendix. We define the corresponding contributions to the zeta function 
Cp(s), Eq. (|1 to be 


An,p{s) — 


sinTTS 


E(^^ +1) 


dz 


1=1 


/21/\-2s 9 

(■r) Yz 


(64) 
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(n = —1,0,1,2,3). The contribution from n = — 1 can be obtained doing the same calculations as in for the 
massive scalar field in a ball with the result 


^-i,p(s) 


(s-i) 
2v^r(s +1) 


Ch{2s 


2; 3/2) 




(65) 


The contribution following for the residuum at s = —1/2 of the zeta function is 

127 

ResA_i^p(s = -1/2) = - cz). (66) 

The z-integral appearing in ^o,p cannot be done by elementary means. However the analytical structure around 
s = —1/2 is easily found. First, for x —> 0 the integrand behaves like 0{z) at s = —1/2 and the integral exists at the 
lower bound. Furthermore one can verify, that for z ^ oo it behaves as 0{z~^) and the integral is well defined also 
at the upper integration limit. Thus Ao,p is analytical around s = —1/2 and gives no contribution to the pole of the 
zeta function. 

Let us now come to the most delicate part of the calculation. As we will see even the pole contribution will be, in 
general, a very involved integral which can only be analysed numerically. This will happen for the contribution arising 
from D^^p{z). But let us continue step by step. First, the angular momentum sum in the Ai^p can be performed in 
terms of the Hurwitz zeta function, 

qi'n TT PCC) O 

A„,p(s) = 2^- -R^‘^CHi2s + n - 1; 3/2) / dzz-^^^D^,p{z) (67) 

TT Jq dZ 

In the limit z ^ 0 all Dn,p{z) converge to a constant so that the integrals always exist at the lower bound. At the 
upper bound however the integral might not be defined at s = —1/2. This happens for example for Di^p{z). Here one 
finds 


-Di.p(z) 


l(ci-C2)l+0 



( 68 ) 


and the integral as it stands does not exist at s = —1/2. However, the analytical continuation is obtained by writing 
for Ai^p instead of (p^ the following. 


Ai,p(s) =2^^R2 s^^^(2s; 3/2) rizz”"* 
^ Jo 


Di,piz) + 


Cl - C 2 1 \ 

8 {z + iyj 


1 s(sin7rs) 

2 sin(27rs) 


i?"*(ci-C2)Cff(2s;3/2). 


(69) 


Here, the asymptotics for z oo has been added and subtracted, the first term is now analytical around s = —1/2, 
the second one clearly contains the pole, namely 

ResAi,p(s =-1/2) =-^^^(ci-C 2 ). (70) 

For A2 ,p similar arguments than the ones provided for Ao,p show that there is no pole. Finally we are left with 

qi'n TT S POO O 

A3.p(s) = 2^- -R^^(:h{2s + 2;3/2) dzz-^^ (71) 

TT Jo 

Here at s = —1/2 the Hurwitz zeta function has already a pole such that the residue of A^ p is determined by the 
integral. In detail, after a partial integration one finds 

ResA3,p(s = -1/2) = ^ y dz D^^p^z). (72) 

An integration of D^^p{z) in terms of elementary functions seems impossible. For the case of /li = 1 simplifications 
occur, but still a complete analytical treatment seems only possible for p = 1, this is for the TE modes. For p = —1 
the integral cannot be expressed in terms of usual special functions. Instead we present the behaviour of the complete 
pole (that is for the pole of the complete zeta function) as a function of the velocities of light in Figure 1. 
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FIG. 1. The dependence of the heat kernel coefficient 02 on the speed of light ci 


1 / inside the ball while C 2 = 1 outside 


In the limit of small differences of the velocities of light (dilute dielectric ball) one might expand the residue in 
powers of Ci — C2. Surprisingly, the leading 0{ci — C2) pole coming from A_i and Ai is exactly cancelled by the 
contribution coming from A^, the second order is zero so that we are left with a pole of the order (ci — 02)^. This 
result is in agreement with where the calculation for a dilute ball was performed up to the order (ci — 02)^ and a 

finite result was obtained in the zeta function scheme. Here we see that in the next order of the dilute approximation 
a pole appears. In detail we find 


ResC(s 


-1/2) 


166 (ci - 02)^ 
50057ri? 4 ’ 


(73) 


especially, for equal velocities of light in the interior and exterior region no pole is present. By means of (H), the heat 
kernel coefficient is 


02 = -167r2ResC(s = -1/2) = + OUci - C 2 )*), (74) 

a value quite close to 02, eq. ( p5| ) for the delta sphere. 

Let us mention that also for the T E and T M modes separately the divergent term is of third order in the difference 
of the velocities of light. 

The peculiarity of this divergent term is, that it has the same 1/i? dependence as the finite part of the result. This 
does not happen to quadratic order where independently of the scheme used [^,^,^, no such divergence appears. 
Thus in contrast to the second order calculation, in the next order one must raise again the question how to fix the 
finite part of the Casimir energy. Until this question is not answered, the calculation of finite parts in higher orders 
of the dilute approximation makes no sense. The situation is even more involved for smooth dielectrics as seen in eq. 
( pO| ) and as discussed further in the Conclusions. 

However, before we come to the summary of our main results let us comment on other possible regularization 
schemes in order to make contact with the results of [Hi- It is known, that the zeta function scheme often gives 
finite answers whereas using other techniques divergences show An example is the dilute ball to second order in 
(ci — C2) where the zeta-function method yields a finite result [ |l3| whereas cutoff functions of different types yield 
poles proportional to the volume and surface for example. These divergent terms are also easily recovered by our 
calculation. There are well established connections between zeta function and cutoff regularization of Casimir energies. 
Using a frequency cutoff A the divergent terms for A ^ 0 are 


p.div ^ 3ao 01/2 Qi 

27r2A4 47r3/2A3 871^2^2 


a2 

167r2 


In A 


(75) 
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with the heat-kernel coefficients Oi of the problem at hand. These are immediately determined by the formulas of 
Ai^p already presented. One gets 


O Q 

Oo = --kR 


01/2 


= 


2^3/2. 


1 

C2 

,212 


Oi =- 'kR 

3 


{cl-cl) 
clcl(cl + cl) 


\Cl C2 
,2 „ 2\2 


SttR 


dy [■Di,i(y) -f £>1 _i(y)] 
y oy 


_^3/2 (Cl C 2 ) 

“3/2-’^ (,2+,2)2 

and 02 is the sum of eqs. (|^, ( [^ and (fz^). It is clearly seen, that volume, surface and more divergencies are 
present. Using an expansion in <5 = ci — C2, one finds explicitly 


C jpdiv I cren 

^0 — + ^0 

9 V ,2( V 

2Tr^cl A'^ VTT^cf 


1 S 
dTTcf A3 


23 1A 

153671^7 


+ 0(S^), 


an expansion in parallel with Barton |^, and where £0®” has been taken from that reference. Let us mention that the 
factors of the divergent pieces can not be compared to 0 because there a wave number cutoff has been used. The fact 
that there is no 1/A term is generic in this cutoff definition |^. In contrast, the reason for 1/A^ not being there is a 
result of cancellations between the TE and TM modes. In the approximation given, volume and surface divergencies 
are absorbed by counterterms and a unique finite result is obtained. However, let us repeat and emphasize that at 
third order a term {S^/R) In A appears and renormalization and normalization of this term remains unclear. 


V. CONCLUSIONS 

We have calculated the lowest heat kernel coefficients for the delta sphere and for the dielectric ball. For a massless 
field, it is the coefficient 02 which determines the interpretation of the ground state energy . In both cases it does 
not vanish. For a delta sphere it is the simple expression (^bj), for the dielectric ball it is an involved expression, 
including the integral ([^. Therefore, just subtracting the “empty space” or “unbounded medium” contribution is 
clearly insufficient to make the ground state energy finite. At this point the conventional field theoretical approach 
fails to give a unique result. It must be supplemented by an extra normalization condition, which is not at hand by 
now. 

A common feature of the two cases considered here is that for small coupling to the background, i.e., for small a 
in or e ^ 1 in (|5^), the lowest two orders vanish ( contrary to the third and higher orders). Thereby nontrivial 
cancellations appear in intermediate expressions. Therefore, the good news is that calculations in the leading orders 
of the dilute dielectric ball approximation are physically relevant and give a unique result. It seems to be true for all 
penetrable boundary conditions . This observation is supported by an example of the paper Q , where Casimir-Polder 
forces between the particles constituting a dielectric sphere, had been summed up over the volume. In fact, the 
potential had been integrated using dimensional regularization. It turned out, that the analytic continuation to the 
physical dimension yields a unique and finite value. In general, one expects the Casimir-Polder forces to be equivalent 
to the ground state energy . Because these forces are nonadditive, the result found in should have relevance to a 
dilute dielectric medium and in fact the identity of the van der Waals force and the Casimir effect has been shown 
in [^. Thus to the leading non-vanishing order the situation is quite clear now for the dilute dielectric ball, but, as 
we have described in detail at the end of section 4 in higher orders the situation remains unclear. Even worse, if 
we consider a smoothly varying dielectric background, we have seen that already in the second order an ambiguity 
proportional to 02, eq. (pO|), remained. If one is going to calculate the Casimir energy by summing up the van der 
Waals energies between the molecules that make up the medium taking into account that e and y may vary in space, 
we expect that in the scheme used in [Q in contrast to the dilute calculation divergencies or at last a logarithmic 
contribution delivering a nonuniquness will show up when the regularization is removed. 

So for the case of nonconstant dielectricum the way how to fix a unique finite value of the Casimir energy still 
remains unclear and further investigations will be necessary to shed more light into this very basic question. 
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APPENDIX A: ASYMPTOTIC FUNCTIONS Dn{Z) 

In this appendix we give a list of the functions used in order to get the analytical properties of the relevant zeta 
function. We have. 


2A{c2ti+Cit2^P)^ 

I—02^1 (—3ti-l-5ti^-l-3t2 + 7t2^)+cit2 (3ti-l-7ti^ — 3^2+ 512^) 


(Al) 


16 (c2 fl -I- Cl t2£,P) 

{c\t\ {t\ — &t\ +ht\ —t\-\-— 1 

-1-4 Cl C 2 1\ ^P + c\ [-t\ + 6 ti-7tl + tl- 6 t 2 +5tl) i^P) , 


~ 5760 (C2 h + Cl t2 ^Pf 
{cl tf (375 tf - 4779 ^ + 9945 tj - 5525 tf 
+tl (345-494Hi-kll655t^-73151^)) 

-3 Cl ci tf t2 (1659 1\ - 3465 tl + 19251? + 120 tl tl - 720 4 tl + 600 1\ tl 
-151? (9 + 8tl-48t^ + 56t^) +4 (-105 + 1581- 3735+ 22751^)) ^p 
+3 cl C2 tl 4 (1581 4 - 3735 4 + 22751? - 120 4 4 + 720 tj 4 - 840 4 4 
+ 151? (-7 + 8t^-48t^ + 40tf) +4 (-135+ 1659t^-3465t^ + 19251^)) 
+cl 4 (-345 4 + 4941 4 - 11655 4 + 7315 tf 
+4 (-375 + 4779 4 - 9945 4 + 5525 4 )) > 


(A2) 
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